{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "##### 收纳"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "rm <- read.csv('datas/model_data.csv', sep=',',encoding=\"UTF-8\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 92,
   "metadata": {},
   "outputs": [],
   "source": [
    "y = rm[43]      # 已对数化\n",
    "x = rm[1:42]\n",
    "\n",
    "model_data <- cbind(y, x)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 114,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "\n",
       "Call:\n",
       "lm(formula = model_data)\n",
       "\n",
       "Residuals:\n",
       "     Min       1Q   Median       3Q      Max \n",
       "-10.5369  -0.0872   0.0357   0.1502   2.7108 \n",
       "\n",
       "Coefficients:\n",
       "                       Estimate Std. Error t value Pr(>|t|)    \n",
       "(Intercept)           9.936e+00  4.087e-02 243.107  < 2e-16 ***\n",
       "个人等级             -2.256e-03  1.078e-03  -2.093 0.036376 *  \n",
       "主播等级              1.818e-02  1.724e-03  10.542  < 2e-16 ***\n",
       "直播时长              1.129e-02  1.210e-03   9.330  < 2e-16 ***\n",
       "距离上次直播时长     -4.395e-05  8.519e-05  -0.516 0.605904    \n",
       "贡献榜总人数          1.712e-02  1.334e-03  12.836  < 2e-16 ***\n",
       "贡献榜总金额          3.175e-09  5.961e-09   0.533 0.594329    \n",
       "贡献榜人均金额        6.381e-07  2.046e-07   3.118 0.001824 ** \n",
       "贡献榜平均用户等级    6.282e-03  1.678e-03   3.744 0.000182 ***\n",
       "广泛贡献用户数        1.357e-03  2.614e-03   0.519 0.603646    \n",
       "金主数量              3.824e-02  3.624e-03  10.552  < 2e-16 ***\n",
       "贵宾数               -1.239e-02  2.402e-03  -5.158 2.54e-07 ***\n",
       "守护者个数            2.145e-03  1.016e-03   2.110 0.034838 *  \n",
       "高等级用户数占比     -6.993e-02  4.599e-02  -1.521 0.128378    \n",
       "本周成长值           -1.492e-08  3.551e-07  -0.042 0.966492    \n",
       "本周排名             -2.282e-07  1.229e-07  -1.857 0.063329 .  \n",
       "视频平均点赞数       -9.013e-04  8.905e-04  -1.012 0.311539    \n",
       "视频平均评论数        5.310e-03  1.783e-03   2.979 0.002898 ** \n",
       "视频平均分享数       -6.855e-04  3.488e-03  -0.196 0.844226    \n",
       "视频平均时长          2.392e-06  1.105e-05   0.217 0.828587    \n",
       "视频平均观看次数      7.373e-06  4.782e-06   1.542 0.123187    \n",
       "视频平均每周更新频率  9.690e-04  4.496e-03   0.216 0.829379    \n",
       "性别_1.0              1.859e-02  1.725e-02   1.078 0.281251    \n",
       "性别_2.0              2.025e-02  2.078e-02   0.975 0.329797    \n",
       "区域_华东            -8.922e-03  2.498e-02  -0.357 0.721016    \n",
       "区域_华中            -2.398e-02  2.653e-02  -0.904 0.366068    \n",
       "区域_华北            -8.547e-03  2.647e-02  -0.323 0.746825    \n",
       "区域_华南            -2.341e-02  2.691e-02  -0.870 0.384341    \n",
       "区域_未知             3.752e-03  2.728e-02   0.138 0.890627    \n",
       "区域_海外            -1.522e-02  4.061e-02  -0.375 0.707867    \n",
       "区域_西北            -1.227e-02  3.540e-02  -0.347 0.728816    \n",
       "区域_西南             2.378e-04  2.646e-02   0.009 0.992828    \n",
       "粉丝规模_10.50       -4.969e-03  1.580e-02  -0.314 0.753154    \n",
       "粉丝规模_1万以上     -1.568e-02  2.732e-02  -0.574 0.565958    \n",
       "粉丝规模_50.600      -4.442e-02  1.951e-02  -2.276 0.022835 *  \n",
       "粉丝规模_600.1万     -3.861e-03  2.201e-02  -0.175 0.860754    \n",
       "活跃时间_黄金时段     4.830e-03  1.137e-02   0.425 0.670864    \n",
       "直播标题主题_2.0      6.736e-03  1.276e-02   0.528 0.597476    \n",
       "直播标题主题_3.0     -1.415e-02  1.235e-02  -1.146 0.251636    \n",
       "蓝光清晰度_高清晰度  -6.900e-03  1.333e-02  -0.518 0.604796    \n",
       "公告主题_直播时间     1.624e-02  1.256e-02   1.293 0.196037    \n",
       "公告主题_粉丝马甲     5.392e-03  1.462e-02   0.369 0.712236    \n",
       "是否点亮徽章_1.0      5.841e-02  1.804e-02   3.237 0.001210 ** \n",
       "---\n",
       "Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1\n",
       "\n",
       "Residual standard error: 0.5197 on 10918 degrees of freedom\n",
       "Multiple R-squared:  0.315,\tAdjusted R-squared:  0.3123 \n",
       "F-statistic: 119.5 on 42 and 10918 DF,  p-value: < 2.2e-16\n"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "result <- lm(model_data)\n",
    "summary(result)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "heading_collapsed": true
   },
   "source": [
    "###### 逐步回归分析"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 115,
   "metadata": {
    "hidden": true,
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Start:  AIC=-14303.16\n",
      "平均小时人气 ~ 个人等级 + 主播等级 + 直播时长 + 距离上次直播时长 + \n",
      "    贡献榜总人数 + 贡献榜总金额 + 贡献榜人均金额 + 贡献榜平均用户等级 + \n",
      "    广泛贡献用户数 + 金主数量 + 贵宾数 + 守护者个数 + 高等级用户数占比 + \n",
      "    本周成长值 + 本周排名 + 视频平均点赞数 + 视频平均评论数 + \n",
      "    视频平均分享数 + 视频平均时长 + 视频平均观看次数 + 视频平均每周更新频率 + \n",
      "    性别_1.0 + 性别_2.0 + 区域_华东 + 区域_华中 + 区域_华北 + \n",
      "    区域_华南 + 区域_未知 + 区域_海外 + 区域_西北 + 区域_西南 + \n",
      "    粉丝规模_10.50 + 粉丝规模_1万以上 + 粉丝规模_50.600 + 粉丝规模_600.1万 + \n",
      "    活跃时间_黄金时段 + 直播标题主题_2.0 + 直播标题主题_3.0 + \n",
      "    蓝光清晰度_高清晰度 + 公告主题_直播时间 + 公告主题_粉丝马甲 + \n",
      "    是否点亮徽章_1.0\n",
      "\n",
      "                       Df Sum of Sq    RSS    AIC\n",
      "- 区域_西南             1     0.000 2949.3 -14305\n",
      "- 本周成长值            1     0.000 2949.3 -14305\n",
      "- 区域_未知             1     0.005 2949.3 -14305\n",
      "- 粉丝规模_600.1万      1     0.008 2949.4 -14305\n",
      "- 视频平均分享数        1     0.010 2949.4 -14305\n",
      "- 视频平均每周更新频率  1     0.013 2949.4 -14305\n",
      "- 视频平均时长          1     0.013 2949.4 -14305\n",
      "- 粉丝规模_10.50        1     0.027 2949.4 -14305\n",
      "- 区域_华北             1     0.028 2949.4 -14305\n",
      "- 区域_西北             1     0.032 2949.4 -14305\n",
      "- 区域_华东             1     0.034 2949.4 -14305\n",
      "- 公告主题_粉丝马甲     1     0.037 2949.4 -14305\n",
      "- 区域_海外             1     0.038 2949.4 -14305\n",
      "- 活跃时间_黄金时段     1     0.049 2949.4 -14305\n",
      "- 距离上次直播时长      1     0.072 2949.4 -14305\n",
      "- 蓝光清晰度_高清晰度   1     0.072 2949.4 -14305\n",
      "- 广泛贡献用户数        1     0.073 2949.4 -14305\n",
      "- 直播标题主题_2.0      1     0.075 2949.4 -14305\n",
      "- 贡献榜总金额          1     0.077 2949.4 -14305\n",
      "- 粉丝规模_1万以上      1     0.089 2949.4 -14305\n",
      "- 区域_华南             1     0.204 2949.6 -14304\n",
      "- 区域_华中             1     0.221 2949.6 -14304\n",
      "- 性别_2.0              1     0.257 2949.6 -14304\n",
      "- 视频平均点赞数        1     0.277 2949.6 -14304\n",
      "- 性别_1.0              1     0.314 2949.7 -14304\n",
      "- 直播标题主题_3.0      1     0.355 2949.7 -14304\n",
      "- 公告主题_直播时间     1     0.452 2949.8 -14304\n",
      "<none>                              2949.3 -14303\n",
      "- 高等级用户数占比      1     0.625 2950.0 -14303\n",
      "- 视频平均观看次数      1     0.642 2950.0 -14303\n",
      "- 本周排名              1     0.932 2950.3 -14302\n",
      "- 个人等级              1     1.183 2950.5 -14301\n",
      "- 守护者个数            1     1.203 2950.6 -14301\n",
      "- 粉丝规模_50.600       1     1.400 2950.8 -14300\n",
      "- 视频平均评论数        1     2.397 2951.7 -14296\n",
      "- 贡献榜人均金额        1     2.627 2952.0 -14295\n",
      "- 是否点亮徽章_1.0      1     2.831 2952.2 -14295\n",
      "- 贡献榜平均用户等级    1     3.786 2953.1 -14291\n",
      "- 贵宾数                1     7.187 2956.5 -14278\n",
      "- 直播时长              1    23.515 2972.9 -14218\n",
      "- 主播等级              1    30.022 2979.4 -14194\n",
      "- 金主数量              1    30.076 2979.4 -14194\n",
      "- 贡献榜总人数          1    44.506 2993.8 -14141\n"
     ]
    },
    {
     "ename": "ERROR",
     "evalue": "Error in eval(predvars, data, env): 找不到对象'平均小时人气'\n",
     "output_type": "error",
     "traceback": [
      "Error in eval(predvars, data, env): 找不到对象'平均小时人气'\nTraceback:\n",
      "1. step(result)",
      "2. eval.parent(fit)",
      "3. eval(expr, p)",
      "4. eval(expr, p)",
      "5. lm(formula = 平均小时人气 ~ 个人等级 + 主播等级 + 直播时长 + \n .     距离上次直播时长 + 贡献榜总人数 + 贡献榜总金额 + 贡献榜人均金额 + \n .     贡献榜平均用户等级 + 广泛贡献用户数 + 金主数量 + 贵宾数 + \n .     守护者个数 + 高等级用户数占比 + 本周成长值 + 本周排名 + 视频平均点赞数 + \n .     视频平均评论数 + 视频平均分享数 + 视频平均时长 + 视频平均观看次数 + \n .     视频平均每周更新频率 + 性别_1.0 + 性别_2.0 + 区域_华东 + \n .     区域_华中 + 区域_华北 + 区域_华南 + 区域_未知 + 区域_海外 + \n .     区域_西北 + 粉丝规模_10.50 + 粉丝规模_1万以上 + 粉丝规模_50.600 + \n .     粉丝规模_600.1万 + 活跃时间_黄金时段 + 直播标题主题_2.0 + \n .     直播标题主题_3.0 + 蓝光清晰度_高清晰度 + 公告主题_直播时间 + \n .     公告主题_粉丝马甲 + 是否点亮徽章_1.0)",
      "6. eval(mf, parent.frame())",
      "7. eval(mf, parent.frame())",
      "8. stats::model.frame(formula = 平均小时人气 ~ 个人等级 + 主播等级 + \n .     直播时长 + 距离上次直播时长 + 贡献榜总人数 + 贡献榜总金额 + \n .     贡献榜人均金额 + 贡献榜平均用户等级 + 广泛贡献用户数 + 金主数量 + \n .     贵宾数 + 守护者个数 + 高等级用户数占比 + 本周成长值 + 本周排名 + \n .     视频平均点赞数 + 视频平均评论数 + 视频平均分享数 + 视频平均时长 + \n .     视频平均观看次数 + 视频平均每周更新频率 + 性别_1.0 + 性别_2.0 + \n .     区域_华东 + 区域_华中 + 区域_华北 + 区域_华南 + 区域_未知 + \n .     区域_海外 + 区域_西北 + 粉丝规模_10.50 + 粉丝规模_1万以上 + \n .     粉丝规模_50.600 + 粉丝规模_600.1万 + 活跃时间_黄金时段 + \n .     直播标题主题_2.0 + 直播标题主题_3.0 + 蓝光清晰度_高清晰度 + \n .     公告主题_直播时间 + 公告主题_粉丝马甲 + 是否点亮徽章_1.0, \n .     drop.unused.levels = TRUE)",
      "9. model.frame.default(formula = 平均小时人气 ~ 个人等级 + 主播等级 + \n .     直播时长 + 距离上次直播时长 + 贡献榜总人数 + 贡献榜总金额 + \n .     贡献榜人均金额 + 贡献榜平均用户等级 + 广泛贡献用户数 + 金主数量 + \n .     贵宾数 + 守护者个数 + 高等级用户数占比 + 本周成长值 + 本周排名 + \n .     视频平均点赞数 + 视频平均评论数 + 视频平均分享数 + 视频平均时长 + \n .     视频平均观看次数 + 视频平均每周更新频率 + 性别_1.0 + 性别_2.0 + \n .     区域_华东 + 区域_华中 + 区域_华北 + 区域_华南 + 区域_未知 + \n .     区域_海外 + 区域_西北 + 粉丝规模_10.50 + 粉丝规模_1万以上 + \n .     粉丝规模_50.600 + 粉丝规模_600.1万 + 活跃时间_黄金时段 + \n .     直播标题主题_2.0 + 直播标题主题_3.0 + 蓝光清晰度_高清晰度 + \n .     公告主题_直播时间 + 公告主题_粉丝马甲 + 是否点亮徽章_1.0, \n .     drop.unused.levels = TRUE)",
      "10. eval(predvars, data, env)",
      "11. eval(predvars, data, env)"
     ]
    }
   ],
   "source": [
    "tstep <- step(result)\n",
    "summary(tstep)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 112,
   "metadata": {
    "hidden": true
   },
   "outputs": [
    {
     "ename": "ERROR",
     "evalue": "Error in eval(expr, envir, enclos): 找不到对象'tstep'\n",
     "output_type": "error",
     "traceback": [
      "Error in eval(expr, envir, enclos): 找不到对象'tstep'\nTraceback:\n"
     ]
    }
   ],
   "source": [
    "tstep"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "hidden": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "hidden": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 98,
   "metadata": {
    "code_folding": [
     0
    ],
    "hidden": true
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Importance of components:\n",
       "                         Comp.1     Comp.2     Comp.3     Comp.4     Comp.5\n",
       "Standard deviation     2.581598 1.68899896 1.57259668 1.41655212 1.24890360\n",
       "Proportion of Variance 0.158682 0.06792184 0.05888239 0.04777666 0.03713715\n",
       "Cumulative Proportion  0.158682 0.22660389 0.28548628 0.33326294 0.37040009\n",
       "                           Comp.6     Comp.7     Comp.8    Comp.9    Comp.10\n",
       "Standard deviation     1.21066774 1.16104631 1.13987840 1.1242881 1.11389067\n",
       "Proportion of Variance 0.03489801 0.03209592 0.03093626 0.0300958 0.02954172\n",
       "Cumulative Proportion  0.40529810 0.43739402 0.46833027 0.4984261 0.52796780\n",
       "                          Comp.11    Comp.12    Comp.13    Comp.14    Comp.15\n",
       "Standard deviation     1.09538453 1.07931164 1.07219578 1.07021412 1.06573684\n",
       "Proportion of Variance 0.02856827 0.02773604 0.02737152 0.02727043 0.02704274\n",
       "Cumulative Proportion  0.55653607 0.58427211 0.61164363 0.63891406 0.66595680\n",
       "                          Comp.16    Comp.17    Comp.18    Comp.19    Comp.20\n",
       "Standard deviation     1.04984337 1.02310110 1.01514497 0.98662379 0.96630062\n",
       "Proportion of Variance 0.02624217 0.02492228 0.02453617 0.02317682 0.02223183\n",
       "Cumulative Proportion  0.69219897 0.71712125 0.74165743 0.76483425 0.78706608\n",
       "                          Comp.21    Comp.22    Comp.23    Comp.24    Comp.25\n",
       "Standard deviation     0.93616635 0.92203341 0.88958271 0.84025643 0.82798422\n",
       "Proportion of Variance 0.02086684 0.02024156 0.01884184 0.01681026 0.01632281\n",
       "Cumulative Proportion  0.80793292 0.82817448 0.84701633 0.86382659 0.88014939\n",
       "                          Comp.26    Comp.27    Comp.28    Comp.29    Comp.30\n",
       "Standard deviation     0.81577229 0.80030940 0.77050652 0.68611170 0.66586659\n",
       "Proportion of Variance 0.01584487 0.01524988 0.01413525 0.01120832 0.01055663\n",
       "Cumulative Proportion  0.89599426 0.91124414 0.92537939 0.93658771 0.94714433\n",
       "                           Comp.31     Comp.32     Comp.33     Comp.34\n",
       "Standard deviation     0.638126289 0.579705975 0.539308692 0.506044871\n",
       "Proportion of Variance 0.009695361 0.008001405 0.006925092 0.006097176\n",
       "Cumulative Proportion  0.956839693 0.964841098 0.971766190 0.977863367\n",
       "                           Comp.35     Comp.36     Comp.37     Comp.38\n",
       "Standard deviation     0.495238313 0.438448110 0.360433622 0.341235801\n",
       "Proportion of Variance 0.005839547 0.004577065 0.003093152 0.002772426\n",
       "Cumulative Proportion  0.983702914 0.988279980 0.991373132 0.994145557\n",
       "                           Comp.39     Comp.40     Comp.41      Comp.42\n",
       "Standard deviation     0.302589087 0.290156127 0.237189896 0.1177998785\n",
       "Proportion of Variance 0.002180004 0.002004538 0.001339501 0.0003304003\n",
       "Cumulative Proportion  0.996325561 0.998330099 0.999669600 1.0000000000\n",
       "\n",
       "Loadings:\n",
       "                     Comp.1 Comp.2 Comp.3 Comp.4 Comp.5 Comp.6 Comp.7 Comp.8\n",
       "个人等级              0.144         0.220                                   \n",
       "主播等级              0.259         0.239                                   \n",
       "直播时长                                  -0.301        -0.347 -0.337 -0.193\n",
       "距离上次直播时长                          -0.193        -0.399 -0.161       \n",
       "贡献榜总人数          0.326                                    -0.124  0.192\n",
       "贡献榜总金额          0.280        -0.211  0.296        -0.147  0.101 -0.172\n",
       "贡献榜人均金额        0.276        -0.183  0.245        -0.110  0.115 -0.207\n",
       "贡献榜平均用户等级    0.213  0.151  0.142 -0.193         0.363        -0.258\n",
       "广泛贡献用户数        0.155                              0.120 -0.132  0.206\n",
       "金主数量              0.335        -0.146  0.143                            \n",
       "贵宾数                0.336                                                 \n",
       "守护者个数            0.268                                    -0.137  0.243\n",
       "高等级用户数占比      0.165  0.125        -0.176         0.398        -0.299\n",
       "本周成长值            0.252        -0.186  0.287        -0.150        -0.184\n",
       "本周排名             -0.118                0.125         0.114  0.222       \n",
       "视频平均点赞数        0.129 -0.514        -0.136                            \n",
       "视频平均评论数        0.155 -0.345                                          \n",
       "视频平均分享数        0.111 -0.480        -0.135                            \n",
       "视频平均时长                        0.203               -0.191              \n",
       "视频平均观看次数      0.120 -0.485        -0.131                            \n",
       "视频平均每周更新频率                0.153               -0.214              \n",
       "性别_1.0                            0.321  0.390               -0.277 -0.100\n",
       "性别_2.0              0.138  0.136 -0.225 -0.391                0.164       \n",
       "区域_华东                           0.112  0.121                       0.251\n",
       "区域_华中                                                                   \n",
       "区域_华北                                                      -0.175 -0.163\n",
       "区域_华南                                                             -0.112\n",
       "区域_未知                          -0.385 -0.242                0.154  0.145\n",
       "区域_海外                                                                   \n",
       "区域_西北                                                                   \n",
       "区域_西南                                                             -0.157\n",
       "粉丝规模_10.50                             0.135         0.101 -0.312  0.115\n",
       "粉丝规模_1万以上      0.189                                    -0.110  0.346\n",
       "粉丝规模_50.600                     0.167               -0.128  0.461 -0.324\n",
       "粉丝规模_600.1万             0.100  0.174 -0.105                            \n",
       "活跃时间_黄金时段                          0.157         0.321  0.313  0.250\n",
       "直播标题主题_2.0                                  0.695                     \n",
       "直播标题主题_3.0                                 -0.693                     \n",
       "蓝光清晰度_高清晰度                                                         \n",
       "公告主题_直播时间    -0.111        -0.341                      -0.181       \n",
       "公告主题_粉丝马甲                   0.267                       0.142  0.149\n",
       "是否点亮徽章_1.0                   -0.185                0.233 -0.208 -0.203\n",
       "                     Comp.9 Comp.10 Comp.11 Comp.12 Comp.13 Comp.14 Comp.15\n",
       "个人等级              0.138  0.211                           0.154   0.189 \n",
       "主播等级                                                                   \n",
       "直播时长                                    -0.119                         \n",
       "距离上次直播时长             0.183           0.128           0.114   0.130 \n",
       "贡献榜总人数                         0.125  -0.126          -0.136         \n",
       "贡献榜总金额                                                               \n",
       "贡献榜人均金额                      -0.111                                 \n",
       "贡献榜平均用户等级    0.123                  0.117                         \n",
       "广泛贡献用户数                       0.161  -0.249  -0.133  -0.238         \n",
       "金主数量                                                                   \n",
       "贵宾数                                                                     \n",
       "守护者个数                           0.148                                 \n",
       "高等级用户数占比      0.166                  0.203   0.142                 \n",
       "本周成长值                          -0.105                                 \n",
       "本周排名                                     0.193                         \n",
       "视频平均点赞数                                                             \n",
       "视频平均评论数                                                             \n",
       "视频平均分享数                                                             \n",
       "视频平均时长                -0.497   0.311   0.161                         \n",
       "视频平均观看次数                                                           \n",
       "视频平均每周更新频率        -0.538   0.311   0.127                         \n",
       "性别_1.0                                                                   \n",
       "性别_2.0                                                                   \n",
       "区域_华东             0.569 -0.305  -0.450  -0.132                         \n",
       "区域_华中            -0.176  0.123   0.274  -0.272   0.717   0.124  -0.241 \n",
       "区域_华北            -0.153  0.151           0.555          -0.537  -0.174 \n",
       "区域_华南            -0.151  0.101   0.121  -0.284  -0.208  -0.230   0.698 \n",
       "区域_未知                   -0.143           0.105                         \n",
       "区域_海外                    0.190                                   0.178 \n",
       "区域_西北                                                                  \n",
       "区域_西南                            0.229          -0.571   0.544  -0.318 \n",
       "粉丝规模_10.50                               0.140                         \n",
       "粉丝规模_1万以上      0.148  0.154   0.213   0.209                         \n",
       "粉丝规模_50.600       0.346  0.105   0.161  -0.170          -0.259  -0.223 \n",
       "粉丝规模_600.1万     -0.423 -0.155  -0.273                   0.184   0.205 \n",
       "活跃时间_黄金时段                    0.114   0.102                         \n",
       "直播标题主题_2.0                                                           \n",
       "直播标题主题_3.0                                                           \n",
       "蓝光清晰度_高清晰度                          0.109  -0.121  -0.108         \n",
       "公告主题_直播时间     0.251          0.257                           0.182 \n",
       "公告主题_粉丝马甲    -0.300         -0.324                          -0.207 \n",
       "是否点亮徽章_1.0            -0.257          -0.321          -0.223  -0.163 \n",
       "                     Comp.16 Comp.17 Comp.18 Comp.19 Comp.20 Comp.21 Comp.22\n",
       "个人等级                              0.169                   0.261   0.167 \n",
       "主播等级                                              0.180   0.167   0.135 \n",
       "直播时长                                                                    \n",
       "距离上次直播时长                             -0.174   0.236  -0.456  -0.336 \n",
       "贡献榜总人数                                                         -0.129 \n",
       "贡献榜总金额                                                                \n",
       "贡献榜人均金额                                                              \n",
       "贡献榜平均用户等级   -0.138          -0.109                                 \n",
       "广泛贡献用户数        0.144                           0.263          -0.444 \n",
       "金主数量                                                                    \n",
       "贵宾数                                                                      \n",
       "守护者个数                                           -0.112                 \n",
       "高等级用户数占比     -0.188          -0.161                  -0.120  -0.128 \n",
       "本周成长值                                                                  \n",
       "本周排名                                      0.103           0.466  -0.640 \n",
       "视频平均点赞数                                                              \n",
       "视频平均评论数                                                              \n",
       "视频平均分享数                                                              \n",
       "视频平均时长                                  0.111                         \n",
       "视频平均观看次数                                                            \n",
       "视频平均每周更新频率                  0.106                                 \n",
       "性别_1.0                                                     -0.158  -0.119 \n",
       "性别_2.0                                              0.209   0.149   0.215 \n",
       "区域_华东             0.127                                                 \n",
       "区域_华中                                    -0.164                         \n",
       "区域_华北             0.231  -0.111                   0.105           0.115 \n",
       "区域_华南            -0.254  -0.182                                         \n",
       "区域_未知                                            -0.272          -0.115 \n",
       "区域_海外                     0.364   0.723   0.209  -0.138  -0.271  -0.124 \n",
       "区域_西北                     0.816  -0.474   0.156   0.141                 \n",
       "区域_西南                    -0.115                                         \n",
       "粉丝规模_10.50       -0.500           0.251  -0.206   0.450   0.150         \n",
       "粉丝规模_1万以上     -0.173          -0.146   0.172  -0.283           0.162 \n",
       "粉丝规模_50.600                                       0.136                 \n",
       "粉丝规模_600.1万      0.524                                                 \n",
       "活跃时间_黄金时段                            -0.136   0.331  -0.495   0.104 \n",
       "直播标题主题_2.0                                                            \n",
       "直播标题主题_3.0                                                            \n",
       "蓝光清晰度_高清晰度           0.247          -0.833  -0.394                 \n",
       "公告主题_直播时间     0.241                                                 \n",
       "公告主题_粉丝马甲    -0.330                          -0.126  -0.140         \n",
       "是否点亮徽章_1.0                      0.158          -0.195                 \n",
       "                     Comp.23 Comp.24 Comp.25 Comp.26 Comp.27 Comp.28 Comp.29\n",
       "个人等级              0.184   0.228   0.262   0.591   0.255   0.198         \n",
       "主播等级                     -0.134                  -0.278   0.159  -0.271 \n",
       "直播时长              0.197  -0.180  -0.136          -0.370   0.414   0.273 \n",
       "距离上次直播时长      0.266                   0.140   0.218  -0.199  -0.187 \n",
       "贡献榜总人数                                                 -0.139         \n",
       "贡献榜总金额                                                                \n",
       "贡献榜人均金额                                                0.131         \n",
       "贡献榜平均用户等级                                                          \n",
       "广泛贡献用户数       -0.290   0.155                   0.307   0.352         \n",
       "金主数量                                                     -0.144         \n",
       "贵宾数                                                       -0.226         \n",
       "守护者个数                   -0.148           0.187  -0.292  -0.396         \n",
       "高等级用户数占比     -0.117                                                 \n",
       "本周成长值                                                    0.125         \n",
       "本周排名              0.386  -0.126                  -0.127                 \n",
       "视频平均点赞数                0.256  -0.116          -0.101                 \n",
       "视频平均评论数               -0.640   0.309           0.279   0.104         \n",
       "视频平均分享数               -0.157                   0.162                 \n",
       "视频平均时长                 -0.180  -0.557   0.263   0.286  -0.117         \n",
       "视频平均观看次数              0.372  -0.176          -0.173                 \n",
       "视频平均每周更新频率          0.239   0.592  -0.227  -0.165                 \n",
       "性别_1.0                                             -0.151   0.129         \n",
       "性别_2.0              0.245                  -0.232   0.115                 \n",
       "区域_华东                                                                   \n",
       "区域_华中                                                                   \n",
       "区域_华北                                                                   \n",
       "区域_华南                                    -0.107                         \n",
       "区域_未知            -0.266           0.118   0.433  -0.167   0.225  -0.161 \n",
       "区域_海外                                    -0.173                         \n",
       "区域_西北                                                                   \n",
       "区域_西南                                                                   \n",
       "粉丝规模_10.50       -0.161                   0.130  -0.103                 \n",
       "粉丝规模_1万以上      0.165   0.193  -0.123  -0.241           0.288  -0.150 \n",
       "粉丝规模_50.600              -0.119           0.131  -0.107                 \n",
       "粉丝规模_600.1万                                     -0.102                 \n",
       "活跃时间_黄金时段     0.273  -0.127                  -0.234   0.287   0.129 \n",
       "直播标题主题_2.0                                                            \n",
       "直播标题主题_3.0                                                            \n",
       "蓝光清晰度_高清晰度                                                         \n",
       "公告主题_直播时间                                                     0.551 \n",
       "公告主题_粉丝马甲                     0.117                           0.566 \n",
       "是否点亮徽章_1.0      0.525                   0.148   0.106          -0.242 \n",
       "                     Comp.30 Comp.31 Comp.32 Comp.33 Comp.34 Comp.35 Comp.36\n",
       "个人等级                                                                    \n",
       "主播等级             -0.224  -0.458  -0.125                           0.225 \n",
       "直播时长              0.124   0.259   0.138                                 \n",
       "距离上次直播时长             -0.152                                         \n",
       "贡献榜总人数                          0.201                          -0.225 \n",
       "贡献榜总金额                                 -0.184                   0.107 \n",
       "贡献榜人均金额                       -0.126  -0.588   0.126  -0.134         \n",
       "贡献榜平均用户等级                                                   -0.655 \n",
       "广泛贡献用户数                       -0.182                           0.111 \n",
       "金主数量                              0.357                           0.154 \n",
       "贵宾数                                0.519                                 \n",
       "守护者个数                    0.221  -0.589  -0.135                         \n",
       "高等级用户数占比                                                      0.582 \n",
       "本周成长值                    0.104  -0.242   0.682  -0.239   0.189  -0.142 \n",
       "本周排名                                                                    \n",
       "视频平均点赞数                                               -0.127         \n",
       "视频平均评论数                                0.141   0.231  -0.385         \n",
       "视频平均分享数                               -0.265  -0.357   0.658         \n",
       "视频平均时长                                                                \n",
       "视频平均观看次数                                      0.131  -0.243         \n",
       "视频平均每周更新频率                                                        \n",
       "性别_1.0                                              0.597   0.361  -0.111 \n",
       "性别_2.0                             -0.102           0.541   0.340         \n",
       "区域_华东                                                                   \n",
       "区域_华中                                                                   \n",
       "区域_华北                                                                   \n",
       "区域_华南                                                                   \n",
       "区域_未知                                             0.135                 \n",
       "区域_海外                                                                   \n",
       "区域_西北                                                                   \n",
       "区域_西南                                                                   \n",
       "粉丝规模_10.50                                                              \n",
       "粉丝规模_1万以上             -0.171                                         \n",
       "粉丝规模_50.600              -0.134                                         \n",
       "粉丝规模_600.1万             -0.144  -0.101                                 \n",
       "活跃时间_黄金时段                                                           \n",
       "直播标题主题_2.0      0.648  -0.269                                         \n",
       "直播标题主题_3.0      0.646  -0.278                                         \n",
       "蓝光清晰度_高清晰度                                                         \n",
       "公告主题_直播时间    -0.146  -0.480                                         \n",
       "公告主题_粉丝马甲            -0.290                                         \n",
       "是否点亮徽章_1.0     -0.104  -0.233  -0.110                                 \n",
       "                     Comp.37 Comp.38 Comp.39 Comp.40 Comp.41 Comp.42\n",
       "个人等级                                                            \n",
       "主播等级                      0.432   0.159                         \n",
       "直播时长                                                            \n",
       "距离上次直播时长                                                    \n",
       "贡献榜总人数         -0.381  -0.183   0.619   0.241                 \n",
       "贡献榜总金额                         -0.234   0.762                 \n",
       "贡献榜人均金额       -0.176           0.175  -0.488                 \n",
       "贡献榜平均用户等级    0.278   0.162                                 \n",
       "广泛贡献用户数        0.122          -0.131                         \n",
       "金主数量              0.730  -0.269   0.150  -0.123                 \n",
       "贵宾数               -0.309   0.263  -0.556  -0.191                 \n",
       "守护者个数                           -0.154                         \n",
       "高等级用户数占比     -0.189  -0.144                                 \n",
       "本周成长值           -0.105                  -0.201                 \n",
       "本周排名                                                            \n",
       "视频平均点赞数                                               -0.751 \n",
       "视频平均评论数                                                      \n",
       "视频平均分享数                                                0.117 \n",
       "视频平均时长                                                        \n",
       "视频平均观看次数                                              0.646 \n",
       "视频平均每周更新频率                                                \n",
       "性别_1.0                                                            \n",
       "性别_2.0                                                            \n",
       "区域_华东                                            -0.458         \n",
       "区域_华中                                            -0.372         \n",
       "区域_华北                                            -0.373         \n",
       "区域_华南                                            -0.357         \n",
       "区域_未知                                            -0.419         \n",
       "区域_海外                                            -0.166         \n",
       "区域_西北                                            -0.199         \n",
       "区域_西南                                            -0.377         \n",
       "粉丝规模_10.50               -0.308  -0.113                         \n",
       "粉丝规模_1万以上             -0.376  -0.185                         \n",
       "粉丝规模_50.600              -0.417  -0.166                         \n",
       "粉丝规模_600.1万             -0.411  -0.181                         \n",
       "活跃时间_黄金时段                                                   \n",
       "直播标题主题_2.0                                                    \n",
       "直播标题主题_3.0                                                    \n",
       "蓝光清晰度_高清晰度                                                 \n",
       "公告主题_直播时间                                                   \n",
       "公告主题_粉丝马甲                                                   \n",
       "是否点亮徽章_1.0                                                    "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# 主成分分析\n",
    "x.pr <- princomp(x, cor=TRUE)\n",
    "summary(x.pr, loadings=TRUE)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "##### AIC准则"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 123,
   "metadata": {},
   "outputs": [],
   "source": [
    "rm <- read.csv('../复赛/datas/modelData.csv', sep=',',encoding=\"UTF-8\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 124,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<style>\n",
       ".list-inline {list-style: none; margin:0; padding: 0}\n",
       ".list-inline>li {display: inline-block}\n",
       ".list-inline>li:not(:last-child)::after {content: \"\\00b7\"; padding: 0 .5ex}\n",
       "</style>\n",
       "<ol class=list-inline><li>'y'</li><li>'粉丝数'</li><li>'个人等级'</li><li>'主播等级'</li><li>'直播时长'</li><li>'贡献榜人均金额'</li><li>'贡献榜平均用户等级'</li><li>'广泛贡献用户数'</li><li>'贵宾数'</li><li>'守护者个数'</li><li>'视频平均评论数'</li><li>'视频平均观看次数'</li><li>'活跃时间_黄金时段'</li><li>'是否点亮徽章_1.0'</li></ol>\n"
      ],
      "text/latex": [
       "\\begin{enumerate*}\n",
       "\\item 'y'\n",
       "\\item '粉丝数'\n",
       "\\item '个人等级'\n",
       "\\item '主播等级'\n",
       "\\item '直播时长'\n",
       "\\item '贡献榜人均金额'\n",
       "\\item '贡献榜平均用户等级'\n",
       "\\item '广泛贡献用户数'\n",
       "\\item '贵宾数'\n",
       "\\item '守护者个数'\n",
       "\\item '视频平均评论数'\n",
       "\\item '视频平均观看次数'\n",
       "\\item '活跃时间\\_黄金时段'\n",
       "\\item '是否点亮徽章\\_1.0'\n",
       "\\end{enumerate*}\n"
      ],
      "text/markdown": [
       "1. 'y'\n",
       "2. '粉丝数'\n",
       "3. '个人等级'\n",
       "4. '主播等级'\n",
       "5. '直播时长'\n",
       "6. '贡献榜人均金额'\n",
       "7. '贡献榜平均用户等级'\n",
       "8. '广泛贡献用户数'\n",
       "9. '贵宾数'\n",
       "10. '守护者个数'\n",
       "11. '视频平均评论数'\n",
       "12. '视频平均观看次数'\n",
       "13. '活跃时间_黄金时段'\n",
       "14. '是否点亮徽章_1.0'\n",
       "\n",
       "\n"
      ],
      "text/plain": [
       " [1] \"y\"                  \"粉丝数\"             \"个人等级\"          \n",
       " [4] \"主播等级\"           \"直播时长\"           \"贡献榜人均金额\"    \n",
       " [7] \"贡献榜平均用户等级\" \"广泛贡献用户数\"     \"贵宾数\"            \n",
       "[10] \"守护者个数\"         \"视频平均评论数\"     \"视频平均观看次数\"  \n",
       "[13] \"活跃时间_黄金时段\"  \"是否点亮徽章_1.0\"  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "colnames(rm)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 125,
   "metadata": {},
   "outputs": [],
   "source": [
    "library(leaps)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 107,
   "metadata": {},
   "outputs": [],
   "source": [
    "# y = rm[43]      # 已对数化\n",
    "# y = exp(y)      # 未对数化\n",
    "# x = rm[1:42]\n",
    "\n",
    "# model_data <- cbind(y, x)\n",
    "\n",
    "# # 重命名列名\n",
    "# colnames(model_data)[1] <- 'y'"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "AIC"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 126,
   "metadata": {},
   "outputs": [],
   "source": [
    "exps <- regsubsets(y~., data=rm, nbest=1, really.big=T)\n",
    "expres <- summary(exps)\n",
    "\n",
    "res_aic <- data.frame(expres$outmat, adjr2=expres$adjr2)\n",
    "# res_aic;       # 查看结果(看不清)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 127,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0gAAANICAMAAADKOT/pAAAAPFBMVEUAAAAaGhozMzNNTU1m\nZmZoaGh8fHyMjIyampqnp6eysrK9vb3Hx8fQ0NDZ2dnh4eHm5ubp6enw8PD////C4q3tAAAA\nCXBIWXMAABJ0AAASdAHeZh94AAAgAElEQVR4nO2djZqrLLaEnfHE3xhHvf97PYImMUm3UDaV\nkE29M8/+evcmUAuoqIjLbBJC/Jns0wKE+BeQkYQIgIwkRABkJCECICMJEQAZSYgAyEhCBEBG\nEiIAMpIQAZCRhAiAjCREAGQkIQIgIwkRABlJiADISEIEQEYSIgAykhABkJGECICMJEQAZCQh\nAiAjCREAGUmIAMhIQgRARhIiADKSEAGQkYQIgIwkRABkJCECICMJEQAZSYgAyEhCBEBGEiIA\nMpIQAZCRhAiAjCREAGQkIQIgIwkRABlJiADISEIEQEYSIgAykhABkJGECICMJEQAZCQhAiAj\nCREAGUmIAMhIQgRARhIiADKSEAGQkYQIgIwkRABkJCECICMJEQAZSYgAyEhCBEBGEiIAMpIQ\nAZCRhAiAjCREAGQkIQIgIwkRABlJiADISEIEQEYSIgAykhABkJGECICMJEQAZCQhAiAjCREA\nGUmIAMhIQgRARhIiADKSEAGQkYQIgIwkRABkJCECICMJEQAZSYgAyEhCBEBGEiIAMpIQAZCR\nhAiAjCREAGQkIQIgIwkRABlJiADISEIEQEYSIgAykhABkJGECICMJEQAZCQhAiAjCREAGUmI\nAMhIQgRARhIiADKSEAGQkYQIgIwkRABkJCECICMJEQAZSYgAyEhCBEBGEiIAMpIQAZCRhAiA\njCREAGQkIQIgIwkRABlJiADISEIEQEYSIgAykhABkJGECICMJEQAZCQhAiAjCREAGUmIAMhI\nQgRARhIiADKSEAGQkYQIgIwkRABkJCECICMJEQAZSYgAyEhCBEBGEiIAMpIQAZCRhAiAjCRE\nAGQkIQIgIwkRABlJiADISEIEQEYSIgAykhABkJGECICMJEQAZCQhAiAjCREAGUmIAMhIQgRA\nRhIiADKSEAH4oJEyIT5MwNkcriq46Wg6iKkEg9PT9CCZSqhRYpXvNhyuKrhpIvEoweD0ND1I\nphJqlFjluw2Hqwpumkg8SjA4PU0PkqmEGiVW+W7D4aqCmyYSjxIMTk/Tg2QqoUaJVb7bcLiq\n4KaJxKMEg9PT9CCZSqhRYpXvNhyuKrhpIvEoweD0ND1IphJqlFjluw2Hqwpumkg8SjA4PU0P\nkqmEGiVW+W7D4aqCmyYSjxIMTk/Tg2QqoUaJVb7bcLiq4KaJxKMEg9PT9CCZSqhRYpXvNhyu\nKrhpIvEoweD0ND1IphJqlFjluw2Hqwpumkg8SjA4PU0PkqmEGiVW+W7D4aqCmyYSjxIMTk/T\ng2QqoUaJVb7bcLiq4KaJxKMEg9PT9CCZSqhRYpXvNhyuKrhpIvEoweD0ND1IphJqlFjluw2H\nqwpumkg8SjA4PU0PkqmEGiVW+W7D4aqCmyYSjxIMTk/Tg2QqoUaJVb7bcLiq4KaJxKMEg9PT\n9CCZSqhRYpXvNhyuKrhpIvEoweD0ND1IphJqlFjluw2Hqwpumkg8SjA4PU0PkqmEGiVW+W7D\n4aqCmyYSjxIMTk/Tg2QqoUaJVb7bcLiq4KaJxKMEg9PT9CCZSqhRYpXvNhyuKrhpIvEoweD0\nND1IphJqlFjluw2Hqwpumkg8SjA4PU0PkqmEGiVW+W7D4aqCmyYSjxIMTk/Tg2QqoUaJVb7b\ncLiq4KaJxKMEg9PT9CCZSqhRYpXvNhyuKrhpIvEoweD0ND1IphJqlFjluw3/7eN1nuX1+Pi7\ny1Ln+NO/PTRNBIuCqQQD7P5YgmQqoUaJVb7b8J8+XVg1p4ffjbmtc8jtv+XD700TwcJgKsHA\nhyCKIJlKqFFile82/JcPX7K8n/o8u2x/WS7yqqye/6yz6vemiWBxMJVg4GMQRZBMJdQoscp3\nG/7Lh+usm/88Z83md+dV3ipyRysYM7GDmEow0BGIJEimEmqUWOW7Df/lw2VmTtz6rLz/asiK\nRV6+Gin/vWkiWBxMJRjwEMQRJFMJNUqs8t2G//Th16NOkQ3LX5v11K756YPrx3igccQCPARx\nBMlUQo0Sq3y34T99+MVITXa+/rU1qw15+/yJ94DGEQt/GY0PBslUQo0Sq3y34T99+NlI9ixv\n/Wtjlf5+QPrakWWCj0EUQTKVUKPEKt9t+E8ffjbSKR+vf23Nqd1YZe0vn/3ekWVyZBQiCJKp\nhBolVvluw3/5cP5kpMqu4i1/PWXmXuz4dJPpoWkiWBxMJRgHBiGGIJlKqFFile82/JcPL6t2\nw23VbitQy99HODYOHw+SqYQaJVb5bsN/+XBjj0CdXZ+zlW0ELkerUcvfEAcH4tNBMpVQo8Qq\n3234Lx/+cWfDehCqM7PPrr6Z7IemiWBxMJVg4GMQRZBMJdQoscp3G/7Tp09WTWErutW0/lTc\n/+2XpolgYTCVYBwYgxiCZCqhRolVvtvwnz697PBeKno20nT/t1+aJoKFwVSCgQ5AJEEylVCj\nxCrfbThcVXDTROJRgsHpaXqQTCXUKLHKdxsOVxXcNJF4lGBwepoeJFMJNUqs8t2Gw1UFN00k\nHiUYnJ6mB8lUQo0Sq3y34XBVwU0TiUcJBqen6UEylVCjxCrfbThcVXDTROJRgsHpaXqQTCXU\nKLHKdxsOVxXcNJF4lGBwepoeJFMJNUqs8t2Gw1UFN00kHiUYnJ6mB8lUQo0Sq3y34XBVwU0T\niUcJBqen6UEylVCjxCrfbThcVXDTROJRgsHpaXqQTCXUKLHKdxsOVxXcNJF4lGBwepoeJFMJ\nNUqs8t2Gw1UFN00kHiUYnJ6mB8lUQo0Sq3y34XBVwU0TiUcJBqen6UEylVCjxCrfbThcVXDT\nROJRgsHpaXqQTCXUKLHKdxsOVxXcNJF4lGBwepoeJFMJNUqs8t2Gw1UFN00kHiUYnJ6mB8lU\nQo0Sq3y34XBVwU0TiUcJBqen6UEylVCjxCrfbThcVXDTROJRgsHpaXqQTCXUKLHKdxsOVxXc\nNJF4lGBwepoeJFMJNUqs8t2Gw1UFN00kHiUYnJ6mB8lUQo0Sq3y34XBVwU0TiUcJBqen6UEy\nlVCjxCrfbThcVXDTROJRgsHpaXqQTCXUKLHKdxsOVxXcNJF4lGBwepoeJFMJNUqs8t2Gw1UF\nN00kHiUYnJ6mB8lUQo0Sq3y34XBVwU0TiUcJBqen6UESZYN8TLiMRFaCwelpepBE2SAfEy4j\nkZVgcHqaHiRRNsjHhMtIZCUYnJ6mB0mUDfIx4TISWQkGp6fpQRJlg3xMuIxEVoLB6Wl6kETZ\nIB8TLiORlWBwepoeJFE2yMeEy0hkJRicnqYHSZQN8jHhMhJZCQanp+lBEmWDfEy4jERWgsHp\naXqQRNkgHxMuI5GVYHB6mh4kUTbIx4TLSGQlGJyepgdJlA3yMeEyElkJBqen6UESZYN8TLiM\nRFaCwelpepBE2SAfEy4jkZVgcHqaHiRRNsjHhMtIZCUYnJ6mB0mUDfIx4TISWQkGp6fpQRJl\ng3xMuIxEVoLB6Wl6kETZIB8TLiORlWBwepoeJFE2yMeEy0hkJRicnqYHSZQN8jHhMhJZCQan\np+lBEmWDfEy4jERWgsHpaXqQRNkgHxMuI5GVYHB6mh4kUTbIx4TLSGQlGJyepgdJlA3yMeEy\nElkJBqen6UESZYN8TLiMRFaCwelpepBE2SAfEy4jkZVgcHqaHiRRNsjHhP+xqjrP8np8/N3F\n1unWGnYs/9JBTCUY8ADEESRRNsjHhP+tqsKqOT38bswfjJT/3jQRLAymEgx8CKIIkigb5GPC\n/1TVJcv7qc+zy/aX5VZe9/hvj00TweJgKsHAdEcTJFE2yMeE/6mqOuvmP89Zs/ndeStvzMud\npolgcTCVYGC6owmSKBvkY8L/VFWZDfOffbZxy5AVG3llNr5+6tY0ESwOphIMTHc0QRJlg3xM\n+J+qWoVs9RTZcP9rn9W7n+aBxhELmO5ogiTKBvmY8MBGarLz5q+vByTKMP65g96lys3hofhs\nkETZIB8THtZI9izv9tc+q/Y/zQONIxYw3dEESZQN8jHhYY10ysfNX5e1iJ1P80DjiAVMdzRB\nEmWDfEz4n6rKn4xUWefc/prvV04Z0kMdxFSCgemOJkiibJCPCQ+wajfcVu0eBT6s5v3UNBEs\nDqYSjAODEEOQRNkgHxP+p6oaewTqbmtzjwLbrN1vmggWB1MJxoFBiCFIomyQjwkPv7PhdmpX\nZv1+00SwOJhKMDDd0QRJlA3yMeF/q+pk1RS2oltN159Oe3djp5hGlqkEA+z/WIIkygb5mPC/\nVTXa3d9LRS9GcsmkDOmhDmIqwcB0RxMkUTbIx4TreSSyEgxOT9ODJMoG+ZhwGYmsBIPT0/Qg\nibJBPiZcRiIrweD0ND1IomyQjwmXkchKMDg9TQ+SKBvkY8JlJLISDE5P04Mkygb5mHAZiawE\ng9PT9CCJskE+JlxGIivB4PQ0PUiibJCPCZeRyEowOD1ND5IoG+RjwmUkshIMTk/TgyTKBvmY\ncBmJrASD09P0IImyQT4mXEYiK8Hg9DQ9SKJskI8Jl5HISjA4PU0Pkigb5GPCZSSyEgxOT9OD\nJMoG+ZhwGYmsBIPT0/QgibJBPiZcRiIrweD0ND1IomyQjwmXkchKMDg9TQ+SKBvkY8JlJLIS\nDE5P04Mkygb5mHAZiawEg9PT9CCJskE+JlxGIivB4PQ0PUiibJCPCZeRyEowOD1ND5IoG+Rj\nwmUkshIMTk/TgyTKBvmYcBmJrASD09P0IImyQT4mXEYiK8Hg9DQ9SKJskI8Jl5HISjA4PU0P\nkigb5GPCZSSyEgxOT9ODJMqmIiOFHVmmEgxOT9ODJMqmIiOFHVmmEgxOT9ODJMqmIiOFHVmm\nEgxOT9ODJMqmIiOFHVmmEgxOT9ODJMqmIiOFHVmmEgxOT9ODJMqmIiOFHVmmEgxOT9ODJMqm\nIiOFHVmmEgxOT9ODJMqmIiOFHVmmEgxOT9ODJMqmIiOFHVmmEgxOT9ODJMqmIiOFHVmmEgxO\nT9ODJMqmIiOFHVmmEgxOT9ODJMqmIiOFHVmmEgxOT9ODJMqmIiOFHVmmEgxOT9ODJMqmIiOF\nHVmmEgxOT9ODJMqmIiOFHVmmEgxOT9ODJMqmIiOFHVmmEgxOT9ODJMqmIiOFHVmmEgxOT9OD\nJMqmIiOFHVmmEgxOT9ODJMqmIiOFHVmmEgxOT9ODJMqmIiOFHVmmEgxOT9ODJMqmIiOFHVmm\nEgxOT9ODJMqmIiOFHVmmEgxOT9ODJMqmIiOFHVmmEgxOT9ODJMqmIiOFHVmmEgxOT9ODJMqm\nIiOFHVmmEgxOT9ODJMqmEo+R6jzL6/Hxdxdbp7uT4xlZphIMeADiCJIom0o0RipsN54efjfm\nD0bKf286mpFlKsHAhyCKIImyqcRipEuW91OfZ5ftL8ttv3aP//bYdDQjy1SCgemOJkiibCqx\nGKnOuvnPc9Zsfnfe9uuYlztNRzOyTCUYmO5ogiTKphKLkcpsmP/ss41bhqzY9GuZja+fujUd\nzcgylWBguqMJkiibSixGWntw25FFNtz/2mf17qdjGVmmEgxMdzRBEmVTiddITXbe/PX1gMQc\nzeMj+y5Vbg4PxWeDJMqmEq2R7Fne7a99Vu1/OpaRZSrBwHRHEyRRNpVojXTKx81fl7WInU/H\nMrJMJRiY7miCJMqmEouR8icjVdY5t7/m+5XHM7JMJRiY7miCJMqmEouRllW74bZq99izD6t5\nPzUdzcgylWAcGIQYgiTKphKLkRp7BOpua3OPPdtm7X7T0YwsUwnGgUGIIUiibCqxGOnHnQ23\nU7sy6/ebjmZkmUowMN3RBEmUTSUWI00n242FrehW0/Wn097d2CmmkWUqwQD7P5YgibKpRGOk\n0e7+Xip6MZKrf+MZWaYSDEx3NEESZVOJxkh/azqakWUqweD0ND1IomwqMlLYkWUqweD0ND1I\nomwqMlLYkWUqweD0ND1IomwqMlLYkWUqweD0ND1IomwqMlLYkWUqweD0ND1IomwqMlLYkWUq\nweD0ND1IomwqMlLYkWUqweD0ND1IomwqMlLYkWUqweD0ND1IomwqMlLYkWUqweD0ND1Iomwq\nMlLYkWUqweD0ND1IomwqMlLYkWUqweD0ND1IomwqMlLYkWUqweD0ND1IomwqMlLYkWUqweD0\nND1IomwqMlLYkWUqweD0ND1IomwqMlLYkWUqweD0ND1IomwqMlLYkWUqweD0ND1IomwqMlLY\nkWUqweD0ND1IomwqMlLYkWUqweD0ND1IomwqMlLYkWUqweD0ND1IomwqMlLYkWUqweD0ND1I\nomwqMlLYkWUqweD0ND1IomwqMlLYkWUqweD0ND1IomwqMlLYkWUqweD0ND1IomwqMlLYkWUq\nweD0ND1Iomwq/4iR/oPAVAIJofKlQRJlU5GRwkKaXgf40iCJsqnISGEhTa8DfGmQRNlUZKSw\nkKbXAb40SKJsKjJSWEjT6wBfGiRRNhUZKSyk6XWALw2SKJuKjBQW0vQ6wJcGSZRNRUYKC2l6\nHeBLgyTKpiIjhYU0vQ7wpUESZVORkcJCml4H+NIgibKpyEhhIU2vA3xpkETZVGSksJCm1wG+\nNEiibCoyUlhI0+sAXxokUTYVGSkspOl1gC8NkiibiowUFtL0OsCXBkmUTUVGCgtpeh3gS4Mk\nyqYiI4WFNL0O8KVBEmVTkZHCQppeB/jSIImyqchIYSFNrwN8aZBE2VRkpLCQptcBvjRIomwq\nMlJYSNPrAF8aJFE2FRkpLKTpdYAvDZIom4qMFBbS9DrAlwZJlE1FRgoLaXod4EuDJMqmIiOF\nhTS9DvClQRJlU5GRwkKaXgf40iCJsqnISGEhTa8DfGmQRNlU4jFSnWd5PT7+7rLW2VdZVg07\nTUczVqTpdYAvDZIom0o0RipsCvXTw+/GfKmzs/+Wjz99bmk6mrEiTa8DfGmQRNlUYjHSJcv7\nqc+zy/aX5fpygnz+t7HM6t+bjmasSNPrAF8aJFE2lViMVGfd/Oc5aza/O69v+ThbC41Z/nvT\n0YwVaXod4EuDJMqmEouRysxcAvVZef/VkBWLkaqsdzUdzViRptcBvjRIomwqsRhpPYnbvmiq\nyIblr6dsavKs+v0SSUb6gS8NkiibSrxGarLz+tcsK+1iw/MnNkQzVqTpdYAvDZIom0q0RrJn\neVcjmcWG6uH66enT0YwVaXod4EuDJMqmEq2RTma1+2okc400PK2NP3w6mrEiTa8DfGmQRNlU\nYjFS/mSkyq7iXY30+G8/NB3NWJGm1wG+NEiibCqxGGlZtRtuq3bb11yXMtIBvjRIomwqsRip\nsUeg7nbTdWuk5d+GrPi96WjGijS9DvClQRJlU4nFSD/ubFgPQvPV0WgWG86/Nx3NWJGm1wG+\nNEiibCqxGGk62cOPPejcz+HWn5r7v/3SdDRjRZpeB/jSIImyqURjpNHu/l4qejbS1BXXf/ul\n6WjGijS9DvClQRJlU4nGSH9rOpqxIk2vA3xpkETZVGSksJCm1wG+NEiibCrvNNJYzRc63Vo4\nqO1kpFe+NEiibCpvNNKY2zWD5U6RjMTmS4MkyqbyRiPVWTu7qc2fVuaCNB3NWJGm1wG+NEii\nbCpvNNK6C2jIT4OMxOdLgyTKpvJGI129MxaFjMTnS4MkyqbyRiOdsuujeadCRqLzpUESZVN5\no5HarFp/uj1EHqzpaMaKNL0O8KVBEmVTeefyd31zT5fJSGy+NEiibCpvvSHb31KbDJWMROZL\ngyTKpqKdDWEhTa8DfGmQRNlU3mykcm/r6R+ajmasSNPrAF8aJFE2lTcbKeyl0b3aaMaKNL0O\n8KVBEmVTebOR7kvgQZGRXvnSIImyqbzZSGNZXNyl8KajGSvS9DrAlwZJlE3l7ad2m1wMAZuO\nZqxI0+sAXxokUTYVGSkspOl1gC8Nkiibipa/w0KaXgf40iCJsqnISGEhTa8DfGmQRNlU3rz7\n+yHz/d7LLMGmoxkr0vQ6wJcGSZRN5ZNG2kuwBTYdzViRptcBvjRIomwqHzy1G8ItOMhIr3xp\nkETZVD55jVSX7jKeTf8XIVSrPwEJAZUw6yYKiUY2lbee2j0QrmEZ6c1RYkKikU1FRgoLc9Yw\n6yYKiUY2lXef2pW5SWx3yStnSajpaMaKOWuYdROFRCObypuNVK8vKO+zoM9TyEhvjRITEo1s\nKp96jEKndrgSZt1EIdHIpvJmI+W3I1LuKIk1Hc1YMWcNs26ikGhkU3n7qV1uHqPo8p1XlB9p\nOpqxYs4aZt1EIdHIpvLuxYZiXbMLdgtpaTqasWLOGmbdRCHRyKby9huy59LYqAvXrG06mrFi\nzhpm3UQh0cim8o/s/o5mrJizhlk3UUg0sqnISGFhzhpm3UQh0cim8ikjXYJeJMlIb40SExKN\nbCrvNlKtLUKHlTDrJgqJRjaVty9/Xwm63CAjvTVKTEg0sqm8/YbseSqyYSiyoFm5ZKS3RokJ\niUY2lQ9sEWrmo1Ef7ulYW200Y8WcNcy6iUKikU3lA0bqzLtkdY2EK2HWTRQSjWwq706iP5/a\nDdlpushIuBJm3UQh0cim8mYjdcZAdptQ0AeSZKS3RokJiUY2lXcvfzemVJWFfRxJRnpvlJiQ\naGRT0c6GsDBnDbNuopBoZFORkcLCnDXMuolCopFNRUYKC3PWMOsmColGNhUZKSzMWcOsmygk\nGtlUZKSwMGcNs26ikGhkU5GRwsKcNcy6iUKikU0lHiPVeZbXT2+YvSx1OveLy0hvjRITEo1s\nKtEYaUnmcHr43ZjbOnsZiV43UUg0sqnEYqRLlvdTnz9uCi+z1UiupwBlpLdGiQmJRjaVWIxU\n2weUzg9Zus7rUah15u6Skd4aJSYkGtlUYjFSmZnX9z0ce4asuBqpdTUdzVgxZw2zbqKQaGRT\nicVIWbb9j6XI1leRlVlXZfne7jwZ6a1RYkKikU0lXiM12Xm6GumnF2U+vCMmmrFizhpm3UQh\n0cimEq2R7Fne8tdsttQ01jsneDLSW6PEhEQjm0q0Rjrl4+OZ3vi0Nv7w6WjGijlrmHUThUQj\nm0osRsqfjFTZVbyHe0c7N5JkpLdGiQmJRjaVWIy0rNoNt1W7H16RKSPFEiUmJBrZVGIxUmOP\nQN3tydmtkfLM7Bwadm7LykhvjRITEo1sKrEY6cedDetBqDb2GuudnJIy0lujxIREI5tKLEaa\nTvcl7vs53PLTmNt/27mRJCO9NUpMSDSyqURjpNHu/l4qejKS/bfT3u4GGemtUWJCopFNJRoj\n/a3paMaKOWuYdROFRCObiowUFuasYdZNFBKNbCoyUliYs4ZZN1FINLKpyEhhYc4aZt1EIdHI\npiIjhYU5a5h1E4VEI5uKjBQW5qxh1k0UEo1sKjJSWJizhlk3UUg0sqnISGFhzhpm3UQh0cim\nIiOFhTlrmHUThUQjm4qMFBbmrGHWTRQSjWwqMlJYmLOGWTdRSDSyqchIYWHOGmbdRCHRyKYi\nI4WFOWuYdROFRCObiowUFuasYdZNFBKNbCoyUliYs4ZZN1FINLKpyEhhYc4aZt1EIdHIpiIj\nhYU5a5h1E4VEI5uKjBQW5qxh1k0UEo1sKjJSWJizhlk3UUg0sqnISGFhzhpm3UQh0cimIiOF\nhTlrmHUThUQjm4qMFBbmrGHWTRQSjWwqMlJYmLOGWTdRSDSyqchIYWHOGmbdRCHRyKbyjxgp\nR2AqgYRQlRDBgkyiS2SkwCQxa7Agk+gSGSkwScwaLMgkukRGCkwSswYLMokukZECk8SswYJM\noktkpMAkMWuwIJPoEhkpMEnMGizIJLpERgpMErMGCzKJLpGRApPErMGCTKJLZKTAJDFrsCCT\n6BIZKTBJzBosyCS6REYKTBKzBgsyiS6RkQKTxKzBgkyiS2SkwCQxa7Agk+gSGSkwScwaLMgk\nukRGCkwSswYLMokukZECk8SswYJMoktkpMAkMWuwIJPoEhkpMEnMGizIJLpERgpMErMGCzKJ\nLpGRApPErMGCTKJLZKTAJDFrsCCT6BIZKTBJzBosyCS6REYKTBKzBgsyiS6RkQKTxKzBgkyi\nS2SkwCQxa7Agk+gSGSkwScwaLMgkuiQiI9V5ltfj4+8u2U8//tB0NGOVxKzBgkyiS+IxUpEZ\nTg+/G/Pshx9/ajqasUpi1mBBJtEl0RjpkuX91OfZZfvLMst++PGnpqMZqyRmDRZkEl0SjZHq\nrJv/PGfN5nfn7OaezY8/Nh3NWCUxa7Agk+iSaIxUZsP8Z5+V918NWXF1z+bHn5uOZqySmDVY\nkEl0STRGWn2ytUuRDde/bn78+dPRjFUSswYLMokuidhITXa+/nXz4/YTG6IZqyRmDRZkEl0S\nr5HsWd7y182Pv306mrFKYtZgQSbRJfEa6ZSP179ufvzt09GMVRKzBgsyiS6Jxkj5k5Equ4pn\n/7r58demoxmrJGYNFmQSXRKNkZZVu+G2ane//Nn++GvT0YxVErMGCzKJLonGSI097HRZfa1M\nRooXLMgkuiQaI/24s2F7PqdTu3jAgkyiS6Ix0nSyx5zCVnSrSUaKEizIJLokHiONdvf3UpGM\nFDdYkEl0STxG+lvT0YxVErMGCzKJLpGRApPErMGCTKJLZKTAJDFrsCCT6BIZKTBJzBosyCS6\nREYKTBKzBgsyiS6RkQKTxKzBgkyiS2SkwCQxa7Agk+gSGSkwScwaLMgkukRGCkwSswYLMoku\nkZECk8SswYJMoktkpMAkMWuwIJPoEhkpMEnMGizIJLpERgpMErMGCzKJLpGRApPErMGCTKJL\nZKTAJDFrsCCT6BIZKTBJzBosyCS6REYKTBKzBgsyiS6RkQKTxKzBgkyiS2SkwCQxa7Agk+gS\nGSkwScwaLMgkukRGCkwSswYLMokukZECk8SswYJMoktkpMAkMWuwIJPokn/FSP+HwFQCCaEq\nIZJEkCAyUliSmGNJBAkiI4UliTmWRJAgMlJYkphjSQQJIiOFJYk5lkSQIDJSWJKYY0kECSIj\nhSWJOZZEkCAyUpt5AxwAACAASURBVFiSmGNJBAkiI4UliTmWRJAgMlJYkphjSQQJIiOFJYk5\nlkSQIDJSWJKYY0kECSIjhSWJOZZEkCAyUliSmGNJBAkiI4UliTmWRJAgMlJYkphjSQQJIiOF\nJYk5lkSQIDJSWJKYY0kECSIjhSWJOZZEkCAyUliSmGNJBAkiI4UliTmWRJAgMlJYkphjSQQJ\nIiOFJYk5lkSQIDJSWJKYY0kECSIjhSWJOZZEkCAyUliSmGNJBAkiI4UliTmWRJAg8RipzrO8\nHh9/d1nrbE+v//bQdDQjm8QcSyJIkGiMVGSG08Pvxnyps7b/lv/uJBnprSQRJEgsRrpkeT/1\neXbZ/rLMbJ19Vs0earPq96ajGdkk5lgSQYLEYqQ66+Y/z1mz+d05W4xULjVnvzcgI72VJIIE\nicVIZTZM5thT3n81ZMWDd2SkWEgiSJBYjJS9HnWKbNj+dcyK3z8dzcgmMceSCBIkXiM12fnB\nV609+dt+YkM0I5vEHEsiSJBojWTP8jZGGvLy9VO3T0czsknMsSSCBInWSCez2n3/65j/fmIn\nI72ZJIIEicVI+ZORKnsidzdScfrhQ/emoxnZJOZYEkGCxGKkZdVuuK3a3S9/zN+GUzHsNh3N\nyCYxx5IIEiQWIzX2CNRl9bWyrZG6nQW7pXQ0I5vEHEsiSJBYjPTjzob11G5w+UhGei9JBAkS\ni5Gmkz38WMfcr4yWn6rtWd7PTUczsknMsSSCBInGSKPd/b1U9GSkTEaKiySCBInGSH9rOpqR\nTWKOJREkiIwUliTmWBJBgshIYUlijiURJIiMFJYk5lgSQYLISGFJYo4lESSIjBSWJOZYEkGC\nyEhhSWKOJREkiIwUliTmWBJBgshIYUlijiURJIiMFJYk5lgSQYLISGFJYo4lESSIjBSWJOZY\nEkGCyEhhSWKOJREkiIwUliTmWBJBgshIYUlijiURJIiMFJYk5lgSQYLISGFJYo4lESSIjBSW\nJOZYEkGCyEhhSWKOJREkiIwUliTmWBJBgshIYUlijiURJIiMFJYk5lgSQYLISGFJYo4lESSI\njBSWJOZYEkGC/CNG+h/C53T+I0C9nUh3y0gCRkZ6RUYSMDLSKzKSgJGRXpGRBIyM9IqMJGBk\npFdkJAEjI70iIwkYGekVGUnAyEivyEgCRkZ6RUYSMDLSKzKSgJGRXpGRBIyM9IqMJGBkpFdk\nJAEjI70iIwkYGekVGUnAyEivyEgCRkZ6RUYSMDLSKzKSgJGRXpGRBIyM9IqMJGBkpFdkJAEj\nI70iIwkYGekVGUnAyEivyEgCRkZ6JR4j1XmW1+P972OVZVVvK77ye9Ma2XciI70SjZEKa5XT\n/Re5/YVx0tVH+e9Na2TfiYz0SixGumR5P/V5drn+os4q80d5K9Hd/+21aY3sO5GRXonFSHXW\nzX+es+b6izwzp3n307kxL3/63Nq0RvadyEivxGKkMhvmP/vsyS3307kyG6dfkZHeioz0SixG\nWg89TwsKddauP/VZvfdpjew7kZFeidlI5+zuntcDUrZFI/tOZKRXYjZSW+bXa6beLD3sfFoj\n+05kpFdiNtJMtZ7bLWsRv39aI/tOZKRXYjFS/rORxnW1Id+vXEZ6KzLSK7EYaVm1G55X7VZn\nvazmPZfSyL4TGemVWIzU2HO37r66sNxHGpa9Du1t9e6XpjWy70RGeiUWI/28s2EsFweVdqvQ\nTtMa2XciI70Si5Gmk13GLmxFtqb8/ov5H3fuxk4y0puRkV6Jxkij3f29VLTUNP/itJ7R7Wz8\nXv5dI/tOZKRXojHS35rWyL4TGekVGUnAyEivyEgCRkZ6RUYSMDLSKzKSgJGRXpGRBIyM9IqM\nJGBkpFdkJAEjI70iIwkYGekVGUnAyEivyEgCRkZ6RUYSMDLSKzKSgJGRXpGRBIyM9IqMJGBk\npFdkJAEjI70iIwkYGekVGUnAyEivyEgCRkZ6RUYSMDLSKzKSgJGRXpGRBIyM9IqMJGBkpFdk\nJAEjI73yjxhJiA8TcDaHqyoQmCKVjlfI15Y+goz0z5WORsjXlj6CjPTPlY5GyNeWPoKM9M+V\njkbI15Y+goz0z5WORsjXlj6CjPTPlY5GyNeWPoKM9M+VjkbI15Y+goz0z5WORsjXlj6CjPTP\nlY5GyNeWPoKM9M+VjkbI15Y+goz0z5WORsjXlj6CjPTPlY5GyNeWPkJ8RhLiC5GRhAiAjCRE\nAGQkIQIgIwkRABlJiADISEIEQEYSIgAykhABkJGECICMJEQAZCTx3YyfFrAQnZEudZFlWVFf\nPMrWv/zs0Qqh9LD8p3WUKh7qu5wc5UElWJd8a+lp6m4/tS4nXZrSJoMsvebUUSIz0vl0S4J5\n6lyFx/z+c5/17trr8/WnfAhdehpWNdX+POjnYvkaYj6H4OMjfyVYl3xr6bn8fd4W1X7J0yax\nauGs+DBRGWkosqLtzTfMeGnmn13zty3Lpl16vXD0UWHmYnHKKvutVOX7X2NY6VX8Ohfycr9c\nuz1oVa45gyoBuuSLS09TOdVLaZft6iw/LyWGLs+w8xaEmIzUZfV2ngx15jgotX3bVvMk66ba\n1fnVqeznCTy2p/lId3Z9sWOlV9GLPTrX9+lipKXK3nmygSmBuuRbSxdN2zbZVCznLOX+AWk+\n+t/Ho8/ynZJ/IyYjlc9ft6Ojk5apO87ng46jgCnalUufd/nJOR+h0sNaYFFzcklp27lkm5dl\nWZ/dp4yYbqhLvrV01c5cTM+ciqF12e4hUX7IrPnP7dBqJtNUw9r7l/rk/s4zRS+5+XLqT6ew\npYdiOeGyalrnd15bV8U09DPn5tSEVIJ1ybeWvp4Ym/907qFM8YgEca7mL/W+bZqqnb+nO+dx\nYP1vX84nBL2r+7HS69mIHVmPk6/tNdI5pBKsS7619NZITXlyFDbXSN0yJulcIxluB9/c/eUx\nH9dvF1Vn1+VUa44BTd3aTu0ca2VY6Xm0zOC31kfOdYnWFF8Zp4vjmAQqQbrkW0vXZT2f2Znz\n47qfLo6z//mSarNqd+LddIrVSIPzdPZibjgMzboyPDi+a9rGTNw8y8uz6cygpU2JxoxsMxc8\nu4qa0/srvfM2CKQE65JvLW06riymdvFb7Vwrv9T2PlJeNqncR+oeXqbmOmjPX2H5dGnHwawQ\nl824/z12/Sofujqfv/uCljbM49laT4yug0aXLRPm+l26Xx5SgnXJt5a2x+axbtfb07yzNYiY\njDRt752dnF8f57ad8rnP6+rUnrb3un8sPF4XhuaSxeDYIoCVtlTXOV85b/V0Qz1PhryrzXfr\nsH9RBepGuuRLS1f2WtQayH4jQTtDeERlpAlYoBzMl/9gurKdL8Kdq1/m5nlhpvhlnPu/ci4J\nYKUN1+uX3uNiucrPbW3XcUunR/2VYF3yraXPpn9zu7dhMI7zOF14B7EZyZtyPXJVZktO7doE\nsaxPn3rb+fO3/P5oYaW7ZrngKdcLn2z/MsncuG+nKlsONY7vU0gJ1iXfWnoy5/xzt5muWG6w\nOcpvSOo+0tn0a+XeaDeZmzbdNFZZMzWuW7eVHZ7OnIPX8/m16+wBKn1dgzut/813z9rrrp4v\nB2b3FT5GgpRgXfKtpSezvLfeQ7AHf6/zhYWUjHRdrXTfDm+adjTzrG1P56ndX7wZmvy6rlb+\nsIXiT6Wv+J2szweiPs/H8+wOcyQr9w9foBKkS762dNffjNTGcVpniM1I5v7ZZG7IZK55eS7n\nL3YzvVr7neacxm05zEPUVfNVydl58wEsvXzEr5i931Qsf9jVhnBKsC751tJV9ojn8JCJzUjX\nHR29c/l7vQfTrP3uvH9jTr8bszOunC9UPb7JsNKTt5HG2nxVmO9SX+d5K8G65FtLb8kieawv\nPiPdzmJ9TmfNl/mw7qAePZ5HakwZU9x9Fw8ubVeq/RjMfLHXPp4fAJRgXfKtpW+4Fz3fRWxG\nqm9HJJ89w0JEQmxGmhp7jXTJiQ8ziuNEcyoVG7EZ6elC0vtzzrOktUDtdz6FlfZp/4nm5YdQ\nSgC664mR/wqyx53vQ4W3Gzw8PubfGUcnFMpXG2n9epzH4OKMY9lMPuajV8RYaXMm6lfuxqn2\n/CCqBOF62LdPybc+1/gXr/QVt/pXS7irHvP72tvFnbMB6O42USNBrPlRzLW71+Ngdk+B19ck\nVnpa77H7M2TdutvS1QCmpFjuCXuqaM22pnoxUuu1NlkiT8ZV1nV92xbuT1Xl9Tb2sH9Dey3u\nr6J/z1XCVxupXb7HLoPPbaR5rOyE9Fp0RkqP5pGhur/hU39dXBPnuBrAdJfLXiXgFrL9cm+n\n8eSju8+AVbL6NNSr+92zrF02zpljnuMbA+7unvg43534jNSV5gBcuk4hetvv82DNP559jFS3\nvbkfXtm55vr2RUqXT6ejHl/tQ7bsV2g9jITpvv3hQzv0yyFxluJ11eF+NP5OXd0zwfgYyTJW\npcsZeHe3Pqna/kp0RiqWM9nMmejHbn3oz+XZPgnmeuS1Xje2tcuWuP3b4Vjp5SObCNylizaz\nszevHMJRJaCR8utWG7/PFN5f7WO/JBpbH473M1J/rhvPSzCsu99BbEZqs2I0Rmqdk3eoi34+\nAZ/Otbkgd21CsRl42tuUOe1e/mKlF933Hz2ep62nZvlI58gmhSoBjdTW69W417lg6V1xk191\nXuxnPIw0ZoXzYLQpfv9RD/b9SD6PqV1b8VlgaZrLfF09VNnZPXn6ahjvE9LV+1jpaRnZYXlc\nyHn6U3e3XZdhlfTn5qrFA/Mch7n6MpccPsWdySzvzKfnZoOtsafd77E/lm01Lnntisy3DaS7\nnyCt3MVmJHta52ukaShNqd7jUmPq5hOZ+ezotJR0fQNjpadlu2Xus2Fl7MzjsWb6LhfMrjkM\nKBmyZenLz0hDfmrN6pfHwxyTfejKmPqWbcK119buKKyK5emS0/6+yeaUVeuBufNccfTu7lcS\nMdJpPSK5N61OZlrZpD12YrqKZvN5e23mrc80w0pPS/t+W/LK5eK4vT6/5AgTUVL1yKldNVt0\nOg+eRqrMIaYt7zgPYkN5s51zF+J4Pl1Xly4nn8sk7+5+JREjrddIHo9RTOZh01tmRmfpZXFq\nOUWqndenWOlpcr8T4crYZOW0ObVzCceUAEayzyv0/dp7zqWvvszRvRXmCsw8Iu9FX67naKPP\neaZ3d7+SiJFuq5s+d9Hm8wfztbTOtX3aZQXMTMjC/fwtVrq79Gf/+xpjeRquRjo7M9Whuidf\nIy0Fl6vLtvZ4T8DF507pA51ZcvBaiLOpZNevicG5awrq7mdSMZK9j5Q5Hhw1xeyodmZ5qO3d\nk2do6+psl8rai3tksdK317R43kUyyfBbm7Jhvrp2TE5QCWykkz3DPC0pLl34L3+v2Bynhcfe\no9HeIVuvec4Od6Dd/UgyRvKkPdl78SYbz6Vx30cqbk/6e0wzrPTD9PWbxH1hF+6aanR8AFWy\nHs9P7mOXrd3s37F7IbbvYNqh9vXoaARc1wI87vSYr8Lm5h9X1tzbH/7fGXdkpGf6er7+reeT\n6+7k7NGLvcZvB78JiZU+MrK9XS67OD+AKrGnOmN/aYvcwxjt9gev5TLfK/zSvM8p9z5sVHYX\nT2eyyY59PzoOvTKSD/N1z/xnXvtcTY61efnDWGSdx32kcu1/vwkJlT4wsvPpy/n+DotgSraf\nczqjr69bfmzdXk94+20jWJ5xvSl2psy8JtU8W9u5hMhIHgz5ehfJb79+Z1ZNB5/7SPMoLYsS\nnjcfkNIP2788lu2n6S7YLRzTfefiPsYMq1i7COZ1sBmAaxLvWX6v09xfd+o40t13EjHSfPlt\nM2jX4KPmPgkTTEoA/+kIlO4f8PvMNTGBh3BM9wYgi4g3gJG8k1jc8eq8I919JxEjQclPhIiF\n2KZrvu6fHLlGUuqBJEEe7wWJzUh1VpizmEvh+TTWgatNg8dOR7uO5F0hUjqaFFIGaM+n932k\ndcYiz6aC96i8k17c3zd//4FAbEa6pSz2G4Mavta0dO4eNU+O+u37gkuvy8L2fcwhU5o43hDz\nI5vsB+5Pd76z5ZzZ9+65M2ncQcoavJNejBv/HJstXkRnpCWJfuF3pKnzvmjaC3zJedp/ZrK5\ntIP5pivbk8d9EKz0tBxF7XVy19a5c2z90zA8PzvqfK+DVX/9wbWRz5wN+z5DN5ym1tzI8jzc\nmZdtgI9DeCe9sPko5jGau9t19/tPxGckf3qzZaWqzBYX3yPYQuWY61Vbn7rW7EJqN9tLw5Se\nbgXWp7BdpzT+aRjs62br2x8X18LnfEg0W+3qi138cjzcZ95GY5T7fWUtG/g8J+78rTlPd7Px\nb12I84jVP+mFTXhRlKfy5PHU2nFiM5J7e/5K05SNfR/I+vfM/xutdh0zusbcoerrycsaWOnp\nVuBsDxjOXe7+2+faxz+cVyhNW3Z2ks1Urgej+pOdkbd8aY6tTeaP0f0W60VHabReN/7Nf7g/\nhiS9uNi72a3X0BwnNiMhb+yzXPvG+7bTULgT5hTmW73p2255/s4xx7DSi+ShapcTJafu40Zy\nn2JuNvE534oxn0ebO9/NEmO/X7xdDi23ZD+7x5i+3Grx2sPhnfTCbFdsUzTSyfv9Amvem7a1\nm4Xd98MXxjr3GSdzumgWA5rlCc+ApZed6mN+auyTvb1zCe+4kdzlN/PWfUCvu7zyPeq2j+9e\nc1ziZaiR/JNemFJdmaCRxrLwXB2uT9YSbVt7P+o/dpXHacO0LMW2eec5baDSlbmemiq76cCZ\ncwxLw/AnIzmXDy/Z1GWDr5Ee/uOiAI2EJL1oq+l6amffTsp6NVlsRtosOjnLjufivGwT61wJ\n14bm0jZ1ffbO9lTnjXn8xmvaQKX7c3Hbldk7l+ygNAzt7SBw8jgM9Of6+gy70eOeYHMXL+/r\n9DRSV3peI5k95f5GwpJe1FX5aCTGtinDFxvJHLnO1/ckOp6UqU2FZeV8B+NKf14e7e7q5WRt\n36ZYafMU2/VQVLjWA6A0DPdXQpceCUqacrm2n7+z565xdUxbmUrLS+VtpO7ksSd/LZ3XlbeR\nsKQXZsn73CR4agdgFo6L+Zu0st2z/9zF/M3Vto1ZJvc6CxzWZ25LuzLcXxwveoVKG0Z7gXe2\nC34u3nGN1LeuhAxn685yOp/8jDTat6m160OyDvq6zovK94gEJr1o2yLFaySEtjYLX8squO8G\nnUud+b1ceb5U6/33T2OlZ//bF/Z5CXnTYkPh/IJZJmI1+MzI1T32hZbugelnEd2p7D2vkcCk\nF/311C4xI3nm/rb/Pn/tXa7j5L0Ptcl9ljPMYlDjefsRLX29w+N1z+xNRhp9MhS19+0TXrcb\n7CB5PMS+XMLk10dqnTfXgaQXi+w6uRuy3rm/X3Kp+3dS73LSeb3MqNbMbPvfqVhpmwBnvb/i\ndWz0T8Nw1EjLYz2N8znWpV6vtKyWay/XHnfWzP/GwntJzTvpxbQ8p29DHFIykn/u76V4V17v\nCwE7NofT/qSpb6tfy/Xsfh4qrHS7WVXzysrjn4ahvi82+FyZGIuaa/wiN3ushtH5hPc6aX2y\nAHRzJ4y595fceT0weidX8U56YefFrUhCRkJyf3dZ25pbrN57tG8Mnq+XGJeTDHdyMKz0bUCx\n94240jDc0gk3ZZn57LMx/78Xc8m+zUj3Ckl/6sbNtHWOz2ykYTmUeh+TvJNeTNM9naR7+8Zh\nYjMSkPu7r9ZpMJbwXbaL33nS1FpvuGcOVvo2+BX4MjmPNAxXhsbvROk+Dz2KL4XdXxSXp6pd\ndLfSvu91udeOHGUSerAPyP1tt6wuXeNKKfiKczKsCQfs0c59LoOVvrV+Rl+lAKVh8Lofek+s\n4DHL1sKee0+wvMLrUHp/JwJJL95BbEaCcn8LEQuxGQnK/Y0Tx7fXlvgU/XNsD3K0r+fojOSb\n+/sY8KmUBbsEw5IPHFPk0/yAJTRzgl1gHE5MMTQekx3q480LAgba62TjMxKTi1/aScNmqJyp\nCh6G1TuxwVJ6P3sEmoahqm+5D7o8cKqPHErICCSmWHZ4LE8tDW1buTqwe0jU0Hu8jDk7la1x\n9gC/UMObtIxUes+t1Q52kFzJZ57yb0AvB3ZkjwDTMNTlZT3E3e+vBaNZHtXzLG3v3folprjm\nu7nulnDMyTEzx/FbR+S7x/T1JXRD31Z5dQ7eJ3diM9Jt1Tv096mhz/xPOKwd1q93V/c30/WR\nUJ+XWT604njqHUzDcLZi50lTNuFPYYA9E/eCPokpnut0zUn7TN9pSQpz6fZ7pD4143KP5FJn\nzBegx2qkgZEg0vmAw4Z1K+fm592y23f2+DfizB7xvOtndwlmrNr20tfO97se5JCRfBJT3P5x\nPC/fFu6q2+sj72O1v1ozdM18bO7ayj5I3fJegR6TkbqHs5jwOcig92S1L3/ulr1uLka2GHtk\nj3g20q7vzuvj3TljsWZNxg4ZyTMxxVzUZDIa1gAct9f7y/WVust2VIeOfv7fpTHfoPYd66zn\n+qIy0nTa+ih4PtISOkGmG8krewS6D3UtOHR14MvqczUCMrDEFO00lsb9We6TdWy5bLz2tfMV\ntfM3iznZvozzed2Qn5NZbKCl/B498zpc4RrJN3vEESOtU6vzyg/pSw/JwBJT3PfkXerccaq2\nlr91t9frQIupPHfnaTD7j1J51JxEZ25O9Lf31bsuI/q8aO8PbreOM5Oy7AdvI2HZI7A0DPYa\n6j5TfB+s96SxXXHxqRRLTPHQYa17C7LpY/u8rmtkJruR1zwfO56G5QXO4D1BfxIxUmWei23L\nO/uH+Lzs200+qf2pcDbpAzJfI2HZI7A0DHah/r7s4fP2S4Bq6bk889oNCySmuK3x24KDc08L\nYqSzNdJpaOqLGcvC4+3tx4jOSM3tQilotX3pykqwwaSaYZ3aYdkj0FO7pflhWZXwv/vsR7fu\nbrhU7hSbE5aYYkvvKr/tbueVtHmIqh3raqzNvaehSMVIzZFVZC8u/pffZprQFxu8skccMVJj\nkmKbr/XAlwPmRuiwTFyvB1eAxBTFdft+V2SlYyGuLpbcR0uwznyi8wHarHecxvmaMR980oof\nIzYjMXd9x7b87ZE9AjLS+WxurZhn3u09lkvodU9zpLje0jy7nQQkpqjMCWNzGcuTPXfd/Qbo\n++3yd+3eUNgvWXjHcRqrtB7so+H/JDNspO1eHv9GnNkjICPZVG9FZvarzdcPjff7CHxom+Uy\nrayXa7X9p4DQxBQmuL7Nr/s2HIfS5XtruUbqHdmPurmu8n4rhfeOt9iMVHrn/j6C906xm4Ue\n33P/e9n+2BYhZ/YIMA2DmWN9VSOHRT+uGV7afHmYfX+HG5qYwqo9l33ttY5hjbQufdTj/iHp\nYr5Ous3XXPDd9iuxGWnIfXN/H8J/s9UyE/2emW53/ubAkT0CTcNwW2xgdKJtfj0CFEETU5ii\nja3RvQ43LXmBbriejrZrHT3xpG4hNiMdOkPyZ/BeutvORNcEfjzPwQ4G7uwRNzzSMLS3VWpD\n4I2ry3FjuRfge5iZvBJTtPeMrB7XMcNjHzszRF+2algkZqRj9+OwtADgYhlU3JWGYVa6OckM\nfJZslY7LNU/YxBQPx0+fe9Xb+2msR/UwYjOSEF/59L2MJL4dr5OX0Tx+2zkTHB8XwapY/Ia9\n3++/6QtbNuA9cLOAXbPD2Y0PsW+k9QrNvG1w6rE8AJAIVsUHeHqs+tNySJgNAlNRmO0HXtMH\nSH4wuV/A/AdM666n7q8s10f+0xYL8ondqXJV3Cy7iQInhNmIYFV8gDSMZFeZWvvG71Neup8D\nsDcrvZIf3POVWi77D3ijjHk++lqjWYt5323wD/IH9qdKe/2PvfRKwkiJ0F5uG1yGS+PcS7tO\nBI/kB4UtXdfrKzHywEen1jza53c/cygmezXi/xDV8h+fDA+vuIxUNI1595JFRvp3aIemvL0e\n3L4seL/48h+P5Ae2uml56Oo8daGf1UeeL2zn5ivcSD4ZHl5xGakyj1Etr66YCCl1VhGsisXP\njOXyUrKqus58B6aAX/KDpbrbfdDQ91fax5dXuMpOXYkZyTPDwysuI617ik+3R54YxGSkl+2N\n7rfIfR91U9s0iG2zTJ/9myZQ8oMnI4XRu62+y68b6ErH/oN1p6KfkbAMD694Gqm17xpjEZOR\nuuzhHVaDM1nVN2LeaFpk5k1z8/ej851dUPKDByP5vaAWYTZSe93SfcmdiYE2/3GAZXh4ZddI\nVXEzEl4zIoJaO8hQZEVrN7aMl6ZwpRX9VtrzZRnboStdlwNQ8oP2mqCgJnzzFvZGzG02Op5e\naLcpL2Z2DzJYhodX9ozUllWCRjJbIu/5uP7Bw5FlvTKy49q7d1H7Jz9oy+XR0XN58tkpDtHU\n/cNig2NhbU2PcT0RdD4aBWR4eGXPSN3m1A6uGBJBrf0Al9qkgi5q5sMUn+U8TkOTn9dx9Ugi\n4Jv8YHtq11dhMzY8p+NyVA6d2lmOZnhwcjNS0/IeNI/QSClgtgi01wnjnmy+yQ8eFxtoy9+j\nRwYk3EhAhgcMu8pQj+ZRwMLvDd+HkJHej9kkdE/FEi75waORgBcb+3E1Uu6RQt2ULTOfNx1d\n8c/wUDftFnfaC/Pa+XUBkZf8W0Z6M2bcm+2ECZf8YDHSOs2m4HfxVyNVnk+Dl6fO84gEZngo\nygdc6zWPB8cmlXRc/zzPbzzaf1ELlPxg3dmwlB7uD50GAtvZcH7a+bdXGMvw8FSp610tT2eZ\nSln8D9GOXX3yT8Tim/zgMj3NsqDrVH1ldq36GmlEBfhneHiq1HXcXaSMzbr1UEb6h1hWvmvf\nK18g+cHjLAs9aWbve7wuwjD8oGYf/yDtTvH7X72/j5azRtZqsIz0AdZJMzZ+owokP3hML0GY\nNNjZIvDKM/8gzcLh7WYjJIeJjCS+DJu7a1mdAPIUspGRxJex2RmCXwQmk/s7EW5rU4PfGvWR\n2x+0tzy6wO3F2QAAB5tJREFUGecLncuaE2wwQXpp8QwSNdJ2ub5xvM36ODLS+6nq2/uGO9fL\n29HkBxsun7uA6LfvaRqn0rGsggXZDpCR+k21TfCHtG7ISG+nLi/rE9td6XqNLJT8YEn3e7t7\n1H9wbNvtJHdlQwczPJj41htPJ5/jS2XuNS2J/4mXVDLS2zHfv+3Ut1XZOAcWSn5gX3Ey3h4C\n/byR+ttfdkEzPHSbdxb4aamWZzo86z+EjPRmxnlQL33teonlFST5wcVaqClXJ33USCa/S2fz\nazXuPddYhocCfMDmenS0f9AuHGWkN3NeH9SZLyF8LsGh5Af5klVhfflWT8uY48a87Wwa7DGj\nLdxvTJ6AIGvUCxsjDflJOxv+KWyyj652v40TSn5QTpfMFK1ttf0HM148XCN1zl08y59+QY5T\nVz1s/3ad3m2MBNwgRpGRPsJ1p03nepweerCnHcp++Xo3Tjp/cP27NaeYK+N0gTZoO6kffLT/\nZMTlvLwpE6j+GDLS+ym2p+qOL1Qs+cFt26rJJsJOA77H/EWxPWK4Hu1FgsSolzNp20q4Wn9A\nRno7w/YmS5bvHzeg5Afj/dGFczdyJ84+F3sbdP7GWLcTBgzyCLfFhpL2jKyM9H6Wk/Zhuefu\nSK2AnfVsblI2nzwgTUtKU2Ok5WWz++ts+IPpIFcj1dgr4xBkpPdjjNSYt1GYdWr3G7w3/3Fh\nThOvF9QfXLMz3Dzt9+bLzX98AJ+iR55JPIiM9GbO9urXLA/bDcwXx/UAmvxguFXYEd/w4qJu\n7AXPctHTeD1hBGV4QG+R1aXvNojDyEhvpjUjWmQ2a3HZNq7rASj5wWgOb12zfl1/MKHZeV09\nsP+9DC4nQUFawCuddQPt4HoD+h+Qkd6PObXrzaYYj5kDJT/Il4dA57o/nxVwVW1N7XqFNBTk\n2N5TpXhuErrToG+68EZGej+3xQaP6Y4kP1jyb5k7/0OZmUT3f9D4V5ZnWKd+SXTsOuwiGR5e\nssdA10sd670uMtL7Mbso7+mk9r9S4eQHzVJh0XfVJ9e/LfbesH1vxf6SyvEMD5M7i9CbkJHe\nz3m8717u/dLoeu9tGdcv/ygSpy+LzWb92+sEDNjAszHSh5f5r8hI4guxexWXdzSFzhV+EBlJ\nfCEmk1Be+iyprAd8+uqLjPQBkLMR+46G3vdULYozukfGri32jxoHMjy0NsmeB8OyuBD6xdSv\nyEgfABhVk3B/KgqzD8LjS/VQbgcGD1oHhywww4PBM0/l9TH2mn9zOpaeT4anGyaXzJXQwJQf\nzZ2YU166LscH2ku7QW6+WI5FjlkGZXiYu/Bi0+773UayAsyGxmJJvq9Mq/8I5ruxnWr7TE3n\ncc7RXup1orXDpckd524fX/FeedLhZSTPDA/brfMed5FMdfbJrDVvg4z0j7Bs5lye8Tx7vA2s\nHZqyuuef8nv1wse56VgCdRkJyvCw3YLqjteUuL/9ltc/MtKbWYx0HVBXnrXRXjA0U1U9PCOx\nW3sMmDNR+/RqZs+nXFvcoQwPsJGG/ulTDGSkN/NkJEfp2jxWZPa3Nssnfd41FgOtFd025Smv\n3bdZsQwPsJE2G5BkpH+GByM5X5lqUjAUWTVO7Xw5ULtKX87xGOnKeHYlwUQzPDzutnOdG7d1\na1cZLjLSv4W95J0HdKx7rx0x7fmyfGMPnestj5e8jsVIzXYDVO3afIBleNjur3Iv27WnzL4M\ntJOR/i1am4ujnc7lKTt5PUdxP/Xpi92LjTaeU7vbptxl87djqQzL8IDRXm8iyUj/FttTu75y\nZGyY7MbPocmv52yF4+wuFiPdWA4YzltDSIYHDHPaWKw/yEj/EI+LDe7lbzuthvZ6uuN8RC4W\nlsu/PLdOGl3HGCTDw8LNm461idtNA/M8YF/7bbY/gIz0Zh6NtD6Mt4fZJHS5pWT9khuyxkLn\na8rXybG9CcvwMNq6bi/GOTve3mx7xDzsePJ4ifwfkJHezGKkZrmynlxfqGbSNNtToy85tbu0\nTV+e7ikfdyc7luFhfZBvPUSfXVsblh7hp1SSkd7MurNhuRQfNq8z+pGX56odM8IY9IImMgjP\nJatmof0pvxrJdWoHZHhYz4aX1RfbhbsslQ30bxgZ6c3csgpfcd9lGbv65Pna4etqWZZ/9D3F\n1/VD146mK1CGh5uRusojyDUlJD0Puoz0AR684xzi5bu3RjNQ1RXv1QserEcCQLVvhoemtnmE\nWiy3Fi3F6oqM9AEeEuc6tyOvthsbdOMybaezD2vbg/9Stm+Gh8HmCD9lefX5pGN3ZCTxbdh1\nhnYazsXp0xeDd2Qk8W3cb632J9/UJ9fTXNqig4z0BdyuvocPZ8bn487wsJ41ro6o/fJf3B4d\nrvTqy1Sp6tvNxy6P5VHysEAZHtZlg8uyWH7xe6XssO5BmlrdkE2UurysNx+70v1AwneCZXjI\nr7kXjDm6qfHLsLmaSUZKFfOF2059W5VNPFfWgcEyPNhHNMzjgjYdg2OxvF6edb9cN/O1rHyS\nMlLUjJXZqVC3xNdxRwCY4cH+sW7Om8/sdrumssXK89VItCe2ZKSoOa8vV51PfjzSJn4rYIYH\n+8cV9+bv9rpL2G5tZKXQ/H+Goyv4+1XFpgAAAABJRU5ErkJggg==",
      "text/plain": [
       "plot without title"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 420,
       "width": 420
      }
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot(exps,scale=\"adjr2\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 128,
   "metadata": {},
   "outputs": [],
   "source": [
    "write.csv(res_aic, file='..\\\\复赛\\\\datas\\\\最优子集_初始版aic.csv')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Cp"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 129,
   "metadata": {},
   "outputs": [],
   "source": [
    "res_cp <- data.frame(expres$outmat, Cp=expres$cp)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 130,
   "metadata": {},
   "outputs": [],
   "source": [
    "write.csv(res_cp, file='..\\\\复赛\\\\datas\\\\最优子集_初始版cp.csv')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "BIC"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 131,
   "metadata": {},
   "outputs": [],
   "source": [
    "res_bic <- data.frame(expres$outmat, BIC=expres$bic)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 132,
   "metadata": {},
   "outputs": [],
   "source": [
    "write.csv(res_bic, file='..\\\\复赛\\\\datas\\\\最优子集_初始版bic.csv')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "##### 其他"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 133,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "\n",
       "Call:\n",
       "lm(formula = rm)\n",
       "\n",
       "Residuals:\n",
       "     Min       1Q   Median       3Q      Max \n",
       "-10.4373  -0.1358  -0.0006   0.1505   3.1364 \n",
       "\n",
       "Coefficients:\n",
       "                     Estimate Std. Error t value Pr(>|t|)    \n",
       "(Intercept)         9.581e+00  5.256e-02 182.280  < 2e-16 ***\n",
       "粉丝数              1.250e-07  2.048e-08   6.106 1.12e-09 ***\n",
       "个人等级           -6.697e-03  1.346e-03  -4.974 6.82e-07 ***\n",
       "主播等级            3.898e-02  2.886e-03  13.510  < 2e-16 ***\n",
       "直播时长            1.216e-02  1.625e-03   7.484 8.75e-14 ***\n",
       "贡献榜人均金额      9.622e-07  1.137e-07   8.463  < 2e-16 ***\n",
       "贡献榜平均用户等级  3.611e-03  1.448e-03   2.494  0.01266 *  \n",
       "广泛贡献用户数      7.767e-03  2.694e-03   2.883  0.00396 ** \n",
       "贵宾数              1.996e-02  1.634e-03  12.213  < 2e-16 ***\n",
       "守护者个数          8.504e-03  8.502e-04  10.002  < 2e-16 ***\n",
       "视频平均评论数      3.990e-03  1.348e-03   2.959  0.00310 ** \n",
       "视频平均观看次数    3.743e-06  8.676e-07   4.314 1.64e-05 ***\n",
       "活跃时间_黄金时段   3.704e-02  1.660e-02   2.231  0.02574 *  \n",
       "是否点亮徽章_1.0    1.217e-01  2.057e-02   5.915 3.58e-09 ***\n",
       "---\n",
       "Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1\n",
       "\n",
       "Residual standard error: 0.4787 on 4229 degrees of freedom\n",
       "Multiple R-squared:  0.4829,\tAdjusted R-squared:  0.4813 \n",
       "F-statistic: 303.8 on 13 and 4229 DF,  p-value: < 2.2e-16\n"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "result <- lm(rm)\n",
    "summary(result)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "R",
   "language": "R",
   "name": "ir"
  },
  "language_info": {
   "codemirror_mode": "r",
   "file_extension": ".r",
   "mimetype": "text/x-r-source",
   "name": "R",
   "pygments_lexer": "r",
   "version": "4.0.5"
  },
  "toc": {
   "base_numbering": 1,
   "nav_menu": {},
   "number_sections": true,
   "sideBar": true,
   "skip_h1_title": false,
   "title_cell": "Table of Contents",
   "title_sidebar": "Contents",
   "toc_cell": false,
   "toc_position": {},
   "toc_section_display": true,
   "toc_window_display": false
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
